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ABSTRACT 

We report on the extremely intense and fast gamma-ray flare above 100 MeV detected by AGILE 
from the Crab Nebula in mid- April 2011. This event is the fourth of a sequence of reported major 
gamma-ray flares produced by the Crab Nebula in the period 2007/mid-2011. These events are 
attributed to strong radiative and plasma instabilities in the inner Crab Nebula, and their properties 
are crucial for theoretical studies of fast and efficient particle acceleration up to 10^^ eV. Here we 
study the very rapid flux and spectral evolution of the event that reached on April 16, 2011 the record- 
high peak flux of F = (26 ± 5) x 10~^ phcm^^ s~^ with a risetime timescale that we determine to be 
in the range 6-10 hrs. The peak flaring gamma-ray spectrum reaches a distinct maximum near 500 
MeV with no substantial emission above 1 GeV. The very rapid risetime and overall evolution of the 
Crab Nebula flare strongly constrain the acceleration mechanisms and challenge MHD models. We 
briefly discuss the theoretical implications of our observations. 



1. INTRODUCTION 

The Crab Nebula is a most remarkable system, consist- 
ing of a rotationally-powered pulsar of large spindown 
luminosity (L^^j ~ 5 x lO'^^ergs"^) interacting with a 
surrounding nebula at the center of the SN1054 super- 
nova remnant (see, e.g., Hester 2008, for a review of 
the Crab properties). The inner nebula is energized by 
the powerful wave/particle output from the pulsar, and 
shows distinctive optical and X-ray brightness enhance- 
ments ("wisps", "knots", and the "anvil" aligned with 
the pulsar "jet") (Scargle 1969; Hester 1995, 2002, 2008; 
Weisskopf 2000) . Despite the small-scale optical and X- 
ray variations detected on timescales of weeks-months, 
the overall high-energy flux resulting from the unpulsed 
synchrotron radiation of the inner Nebula has been con- 
sidered essentially stable for many decades. 

The discovery by the AGILE satellite of a strong 
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gamma-ray flare above 100 MeV from the Crab Nebula in 
September 2010 (jTavani et al. 20101 ITavani et al. 201 1|) 
and the confir ma tion by the Fermi-LAT 
(|Buehler et al. 2010: lAbdo et al."20TT|) started a 
new era of investigation of the Crab Nebula and of the 
particle acceleration processes in general. Three intense 
gamma-ray flaring episodes from the Crab Nebula have 
been report ecp^ in the gamma-ray energy range 100 
MeV - a few GeV by AGILE and Fermi-LAT prior to 
April 2011 (Tavani et al. 2011, hereafter Til; Abdo 
et al. 2011, hereafter All). This activity has been 
attributed to transient emission in the inner Nebula 
due to the lack of any variation in the pulsed (radio 
and high-energy) signal of the Crab pulsar or of any 
detectable alternative counterpart (e.g., Heinke 2010, 
Cusumano et al. 2011) . 

Starting on 2011 April 11-12, a new gamma-ray flar- 
ing episode with substantial emiss ion above 100 MeV 
was first detect ed by Fermi-LAT (jBuehler et al. 20Tll 
[Hays et al. 201l|) a nd then confirmed by AGILE 
(|Tavani et al. 20lT|) . The flare developed in the fol- 
lowing days with substantial gamma-ray emission 2-3 
times the normal average valuep^ until it reached on 
April 16, 2011 the unprecedented high value of = 
(19.6 ±3.7) X 10~^phcm ~^s~^ for a 24-hour integra- 
tion f Striani et al. 20rT|) (hereafter is the Crab pul- 
sar plus nebula flux above 100 MeV). This detection of 

1® Enhanced TeV emission from the Crab nebula has also been 
reported by the air-shower ARGO-YBJ experiment in coincidence 
with the September 2010 event (Aielli et al. 2010). However, 
this claim was not supported by the simult aneous observations by 
VERITAS ( |Ong et al. 20lol l and MAGIC (IMariotti et al. 201011 . 

The Crab steady state (pulsar plus nebula) flux above 
100 MeV detected by AGILE is F^^steady = (2.2 ± 0.1) X 
10'^ phcm~2 s-i. 
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very rapid variations of the gamma-ray emission in the 
April 2011 flare confirms a trend already noticed in the 
Fermi-LAT data for the September 2010 event by Balbo 
et al. 2011. A re-analysis of the September 2010 event 
AGILE data, which will be presented elsewhere, shows 
that the gamma-ray variability on a timescale below 1 
day is confirmed also by AGILE for that event. 

The goal of our paper is threefold: (1) investigate the 
short timescale structure of the gamma-ray emission of 
the April, 2011 event; (2) present the gamma-ray spec- 
trum at the flare peak; (3) briefly discuss the theoreti- 
cal implications of the very rapid variability and overall 
emission. 

2. THE APRIL 2011 GAMMA-RAY SUPER-FLARE 



The AGILE satellite ()Tavani et al. 2009[) has been 
monitoring in spinning moderM the Crab Nebula region 
with optimal exposure during the period February- April, 
2011. Starting on April 10-11, 2011 a noticeable ris- 
ing gamma-ray flux from a source positionally consistent 
with the Cr ab Nebulap^ was reco rded by the Fermi-LAT 
instrument (|Buehler et al. 201ip . Enhanced emission 
with respect to the average flux of the steady pulsar plus 
nebula emission was detected at even larger flux values 
by AGILE above 100 MeV during the following days, 
reaching a value = (6.5 ib 1.5) x 10~^ phcm~^ s~^ 
on April 11-13, 2011 (jTavani et al. 201 ip . Remark- 
ably, the gamma-ray flux increased even more in the 
following days, reaching the 1-day averaged flux of 
= (12.1 ± 0.6) X 10~^phcm-^s~^ on April 14, 
2011 dHays et al. 2011D , and the even larger and cur- 
rently record-breaking 1-day averaged flux of F^ = 
(19.6 ± 3.7) X IQ- ^phcm-'^s-^ on April 15-16, 2011 
(jStriani et al. 20111) . Fig. [l] shows the AGILE 12 hr- 
binned gamma-ray lightcurvcrj (Crab pulsar plus neb- 
ula) during the period 10-19 Aj)ril, 2011 as obtained by 
the standard AGILE maximum likelihood procedure. 

The 12-hr integration between MJD=55667.0 and 
MJD=55667.5 (April 16, 2011) yields for the peak flux 
the value F-y_p — (26 ± 5) x 10~^ phcm~^ s~-^, with a 
pre-trial statistical significance ap = 10. Considering 
the accumulated AGILE exposure on the Crab Nebula in 
spinning mode (equivalent to 850 maps of 12-hr integra- 
tions) , the post-trial significance for the super flare turns 
out to be (Tp^post-t — 9.2. The statistical significance of 
this 12-hr excess emission compared to the steady Crab 
Nebula flux is 7.5 cr. This significance is calculated by 
the likelihood method taking into account the standard 
average flux and position of the Crab pulsar and Nebula, 
together with the catalogued gamma-ray sources in the 
field. After taking into account the accumulated expo- 
sure, we find that the post-trial significance of the excess 
emission is 6.4 a. For a few days in April, 2011, the 

The AGILE spinning mode allows a daily exposure of about 
70 % of the sky depending on solar panel constraints, for an in- 
strument boresight rotation period of about 7 minutes. 

The possibility of a chance positional coincidence with a back- 
ground transient gamma-ray source has been considered in Til and 
All. In light of the results presented also in this paper, we consider 
this probability as negligible. 

We note that a loss of AGILE-GRID telemetry data occurred 
on April 14, 2011 caused by pre-determined satellite tracking activ- 
ity of the AGILE ground station in Malindi (Kenya). This interval 
partially overlaps with the first strong gamma-ray flare detected 
on April 14, 2011. 
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Fig. 1. — Crab (pulsar plus Nebula) gamma-ray 12-hr binned 
lightcurve above 100 MeV detected by AGILE-GRID during the 
period 10-19 April, 2011. The gray horizontal band indicate the 
Crab pulsar plus nebula average flux in the AGILE bandpass, and 
the gray vertical lines mark the time period during which a loss of 
AGILE-GRID telemetry data occurred because of ground station 
activity. Peak emission occurred on April 16, 2011. 



Crab Nebula became the brightest gamma-ray source in 
the sky, rivaling in intensity not only the Vela pulsar, but 
also the brightest blazars ever detected in gamma-rays. 

In order to determine the Crab super-fiare peak time 
and duration, we proceeded as follows. Starting from the 
12- hr binned lightcurve showed in Fig. 1, we obtained 
different lightcurves by progressively shifting their bin 
starting times by 1 hour. These 12 lightcurvef]^ have 
been compared with those calculated from a flux model 
obtained by integrating for each bin the model function 



B + Ae-^l*-*^!/'"'-)'' for t < tp 



where B — 1 y. 10^'"' phcm^^ s^^ day^^ is the flux base- 
line (i.e., the average of the 12 light curves baselines ex- 
cluding the peak emission), tp is the peak time, Gr and 
<Td the rise and decay time constants, respectively, and v 
a measure of the pulse sharpness {v = 1,2 for two-sided 
exponential and gaussian fits, respectively; see, e.g., Nor- 
ris et al. 1996). In general, lower values of v imply a 
more peaked pulse. We also notice that the rise and de- 
cay times, half to maximum amplitude, are obtained as 



Tr.d 



[ln(2)] 



In the following, we conservatively assume a Gaus- 
sian model for the flaring with v — 2, and a^. — ad- 
From the values of a^^d obtained by the best fit of each 
lightcurve (for reduced values ranging from 1.2 and 
1.4), we derived that the super-flare duration ct^ -I- Cd 
is constrained in the 68% confldence level range be- 
tween 14 hr and 26 hr. Taking into account the relation 
Tr = [ln(2)]^/^ (Tr — 0.8 (7^, wc obtaiu the risetime range 

These datapoints were obtained by the AGILE maximum like- 
lihood procedure. 
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Fig. 2. — The AGILE-GRID gamma-ray pulsar-subtracted spec- 
trum of the Crab Nebula super-flare on April 15-16, 2011. The 
AGILE flaring spectral data, marked in red, obtained for a 1-day 
integration (MJD = 245566.4 - 245567.4); data points marked 
in black show the average nebular spectrum (Meyer et al. 2010). 
Pulsar gamma-ray spectral data have been subtracted based on the 
AGILE results presented in Pellizzoni et al. (2009). The red curve 
is the result of the theoretical modelling of the super-flare as dis- 
cussed in the text. The spectral region marked in green shows the 
X-ray spectrum of "source A" which is the most dominant source 
in the Chandra image of the "anvil" region in the inner Crab Neb- 
ula as reported in Til. This flux level is indicative of an X-ray 
upper bound expected from the flare. 

6hr < Tr < 10 hr. This result is independent of the 
lightcurve bin zero-phase {t start) choice. The value of 
the peak time is determined as MJD=55667.3 ± 0.3. 

A complete spectral evolution of the April 2011 event 
will be reported elsewhere. Here we focus on the 1-day 
integrated April 15-16, 2011 super-flare spectrum that 
we show in Fig. [2j Very intense and relatively hard emis- 
sion is detected in the energy range 100 MeV - 1 GeV. 
The AGILE optimal spectral sensitivity in the 50 MeV - 
a few GeV energy range is important in constraining the 
spectrum at relatively low gamma-ray energies. Indeed, 
the 50-100 MeV flux is well constrained by our 95% con- 
fidence level upper limit. The super-flare emission shows 
a very prominent peak of the vF^, spectrum at photon 
energies E^^p ~ 500 MeV. No significant emission is de- 
tected above 1 GeV. 

Fig. [2| shows also the results of our theoretical mod- 
elling of the emission (red curve) that we discuss below. 
Remarkably, the peak power emitted in the gamma-ray 
energy range between 100 MeV and a few hundreds of 
MeV equals the average power emitted in the hard X- 
ray/MeV energy range by the Crab Nebula. 

3. THEORETICAL CONSTRAINTS 

In modelling the April 2011 Crab Nebula gamma-ray 
event we assume that fast and very efficient acceleration 
is occurring at a site in the inner nebula, following the 
discussion of Til and Vittorini et al., 2011 (hereafter 
VI 1). A fraction of the total electron-positron high- 
energy component in the Nebula is impulsively acceler- 
ated at a site of size L. For simplicity, we ignore here 



substantial enhancements due to Doppler boosting, and 
take the Doppler factor S = r^^(l — /3cos6)^^ to be of 
order of a few. 

The physical quantities are constrained within a global 
comparison of a multi-parameter model matching spec- 
tral and timing datgp^ We considered several models 
with the assumption of S in the range 1 — 4 as deduced 
from observations of the South East jet and wisp regions 
(e.g. Hester 2008) . We present here the cases with 6 = 1 
and 5 = 4 as examples of a class of models applied to the 
super-flare spectrum shown in Fig. 3. 

The acceleration process produces, within a timescale 
shorter than any other relevant timescale, a differential 
particle energy distribution (that we model for illustra- 
tion purposes in its simplest form as a single power-law 
distribution dn/d'j = K / / where n is the local 
particle number density, 7 is the particle Lorentz factor 
ranging from ^^in = 10^ to ^raax = 7 x 10^, s = 2 is 
the power-law index, 7b = 5 x 10^, and K is the nor- 
malization factor K = A y. 10^^ cm^-^. For 5 = 1, the 
emitting region has size L = 10^^ cm, and an enhanced 
local magnetic field S/oc = 2 x 10^"^ G that we keep 
constant in our calculations. The total particle num- 
ber required to explain the flaring episode turns out to 
be ^e_/e+ ^ J dV [dn/d-y) ^7 ~ 7 x 10^^, where V is an 
assumed spherical volume of radius L. For 5 = 4, some 
of the physical parameters are slightly different, e.g., 
7,„„^ = 5xl0^ Bloc = l.SxlO-^G, K = 3x10"^° cm'^, 
L = 4 X lO^^^ cm, and iVe^/e+ = 3 x lO'^^ Obviously, the 
physical parameters can differ from these (and are even 
more extreme) for time variations faster than the 1-day 
spectral average of Fig. 2. A discussion of the complete 
lightcurve and physical implications will appear else- 
where. The complex gamma-ray lightcurve shows that 
the acceleration process occurs on a week timescale 
with a succession of short timescale flares, each of which 
has physical parameters which differ from those of the 
peak emission by a factor of a few. 

We flnd that the synchrotron peak photon energy dur- 
ing the flare maximum is Ep^ak = l^^^^^^lliax - 
500 MeV, a value that challenges models of diffusive par- 
ticle acceleration limited by synchrotron cooling, a fact 
already noticed in Til, All, and Vll. The Aprfl 2011 
event confirms even more the extremely short timescale 
of acceleration occurring in the inner Crab Nebula, and 
the existence of a strongly enhanced local magnetic field. 

4. DISCUSSION AND CONCLUSIONS 

The April 2011 Crab Nebula super-flare dramatically 
shows the efficiency of the particle acceleration mech- 
anism operating in the inner nebula. The detected 
gamma-ray luminosity at the peak of the April 16, 2011 
super-flare corresponds to 0.3% of the Crab pulsar spin- 
down luminosity. Despite recent high-resolution observa- 
tions of the inner nebula (see especially the very interest- 

■^^ Having determined from the overall spectral shape the val- 
ues of 76 and index s (see their definitions in the text), we have 
five remaining parameters: 'ymax, the local magnetic field B, the 
electron density A^e, the dimension of the emitting region L, and 
the Doppler factor 5 (the parameter K is derived from A^^ and 
L). These parameters are obtained from the following quanti- 
ties (in the observer frame): the position of the peak emission, 
Ep oc S-y^^^B, the peak emission uF oc S"^ Ne L'^ B^ ■y'^^^, the 
rise time Tr = L/{c5), and the cooling time Tc oc 1/ {B'^^^'yrnax 
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ing sequence of Chandra pointings reported by Tennant 
et al. 2011), there is currently no identification of the 
acceleration site. The anvil region has been suggested 
as a candidate for the September 2010 event (Til), and 
this site may well be active also in the case of the April 
2011 event. More observations and extended monitoring 
of the Crab Nebula is necessary also in light of the " sec- 
ular" X-ray variations as determ ined over a timescale of 
years ( [Wilson- Hodge et al. 2010D . 

The challenge to particle acceleration models applies 
to the main theoretical frameworks that have been so 
far proposed for the Crab Nebula: the magnetohydro- 
dynamical approach of Kennel & Coroniti (1984), dif- 
fusive shock acceleration mechanisms (e.g., Drury 1983, 
Blandford & Eichler 1987), shock drift acceleration (e.g.. 
Kirk 2000, Reville & Kirk 2010), magnetohydrodynam- 
ical instabilities (e.g., Komissarov & Lyubarsky 2004, 



Del Zanna et al. 2004, Camus et al. 2009, Komissarov 
& Lyutikov 2010), ion-driven acceleration (e.g., Arons 
2008). The data presented here contribute in a substan- 
tial way to a step further in deepening our understanding 
of acceleration processes, and may force the previously 
proposed models to be substantially revised. Of partic- 
ular relevance is the possible role of impulsive particle 
acceleration in magnetic field reconnection, and/or run- 
away particle acceleration by transient electric fields vio- 
lating the condition E/B < 1 (that is typically assumed 
in standard models, e.g., deJager et al. 1996). The 
super-flare spectrum of Fig. 2, showing significant emis- 
sion above 200 MeV, implies E/B > 2 (jTavani 2011|) . 
The applicability of these concepts to the Crab Nebula 
flaring activity remains to be tested by future investiga- 
tions. 

Research partially supported by the ASI grant no. 
1/042/10/0. 
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